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; ABSTRACT 


__~ A survey of experimentally known facts is given. The electronic structure of the H,O+ ion 
_has been studied by the LCAO MO method. Three different values of the H-O-H angle have 
_ been considered and the question of the equilibrium angle, which is not known experimentally, 
is discussed. The results give a slight preference to a value of the H—O-H angle just in the 
neighbourhood of 120°, i.e. a very flat molecule ion. The proton affinity of the water mole- 
cule is obtained in satisfactory agreement with the experimental value. Population analysis and 
transformation to two different sets of equivalent orbitals have been made. 


; I. Introduction 


Very few experimental data are known about the structure of the hydronium 
ion H,0*. Most of these have been obtained by investigation of the properties 
of crystals which contain H,0* ions embedded in the crystal lattice. 

The ion itself has, however, for a long time been expected to be a stable 
entity for thermodynamical reasons [1]. By studing the rapid proton transfer in 
water it has been found that the lifetime of the H,0* complex is probably 
about 10~™ sec [2,3] before leaving a proton for an adjacent H,O molecule. 

This time is greater than the characteristic time of the molecular vibrations 
by about a factor of 10. 

Some difficulties in the explanation of the abnormal proton mobility and the 
UV absorption in acidic solutions can be accounted for, according to Wicke, 
Eigen and Ackermann [3,4], by supposing H,0* ions to be surrounded by a 
hydration shell, thus forming a stable pyramidal H,O, complex. 

In 1957 Falk and Giguére [5] gave more direct experimental evidence for the 
existence of the ion in solutions of strong acids. They detected it with infrared 
spectroscopy. Their data support the idea that it exists long enough to undergo 
normal vibrations. 

In the experimental studies discussed below data were not obtained for free 
ions but for ions in some way linked to other molecules. The binding forces, of 
course, also affect the structure of the H,O* to some extent. This influence will 
be discussed later in the paper. It is, however, very interesting for our purposes 
to study the experimental data because they give a good idea of the limits 
within which the molecular parameters may vary. 
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IR and Raman spectra : 

Bethell and Sheppard [6] identified the vibration spectrum of H;0*NO; and 
H,0*Cl0;. Ferriso and Hornig [7] found, in a study of the ionic crystals OH,F, 
OH,Cl, OH,Br, OH;I and OD,Cl the four fundamental frequencies expected for 
a symmetrical pyramid. For their analysis they assumed tetrahedral angles and 
an O-H distance of 0.98 A. 

The same O-H distance was used in an analysis of the vibration spectrum of 
the hydronium ion in hydronium perchlorate, obtained by Taylor and Vidale [8] 
using three different H-O-H angles, 107 °, 111° and 115°. They found that there 
might perhaps be a slight perference for 107 oe ; 


NMR spectra 


With the aid of NMR spectra of some acid monohydrates below — 130°C, — 
Richards and Smith [9] were able to identify the H,O* ion and calculate the — 
O-H distance to 1.02 A, the H-H distance to 1.72 A and the angle H-O-H to 
115°. There is good agreement between the results obtained by them for different 
monohydrates. 

The results of similar studies of the NMR spectrum of hydrogen perchlorate 
monohydrate performed by Kakiuchi et al. [10] gave the results 0.98 A, 1.58 A 
and 110°. ; 


X-ray spectra 


In 1924 an X-ray study of HClO,-H,O was made by Volmer [11]. By com- 
parison with the corresponding X-ray data of ammonium perchlorate he found that _ 
the configuration of the atoms in the two associations were similar, which indicates 
a symmetric top configuration for H,O~ like that of NH,. In later years direct 
observation of the diffraction pattern of the hydrogen atoms in crystals have 
been possible. The effect is very weak, however, and the positions of the H 
atoms can only be determined with a standard deviation of about 0.07 A. Luzzati 
[12, 13] investigated the diffraction pattern in crystalline HNO,-3H,O. He ob- 
tained different individual values for the three O-H distances, depending on 
their positions in the crystal lattice: 1.00, 0.85, and 1.10 A. The angles between 
the O-H:--O hydrogen bonds in which the hydrogen atoms take part are 113°, 
118°, and 103°. The average distance from the hydrogen atoms to the O-O line 
is found to be 0.12 A. 

Young Ku Yoon and Carpenter [14] determined the crystal structure of 
hydrogen chloride monohydrate by X-ray diffraction. They obtained the value 
0.96 +0.08 A for the O-H distance and approximately 117° for the H-O-H angle. 
The distance between a hydrogen atom and the line joining the oxygen atom to 
the corresponding chlorine atom was found to be 0.10 A, which is almost the 
same as in the investigations by Luzzati. 


Theoretical studies 


In 1925 Hund [15] made an attempt to determine the structure of the hy- 
dronium ion by theoretical methods. He estimated with the aid of spectra of 
certain molecules the potential of some atomic ions, which made it possible to 
calculate the probable configuration of H,O and H,O* among others. He found 
for H,O the values 1.03 and 1.09 A for the O-H and H-H distances respec- 
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tively. Corresponding values for H,0* were 1.05 and 1.11 A, which means much 


more peaked forms for the molecules than they have later been believed to have. 
Another result of his calculations was the value 180 (+40) kcal/mole for the 


proton affinity of water, which is in excellent agreement with values found later 


_ (see part V). This agreement possibly has its root in the semi-empirical nature 
of his investigation, where known experimental data are used at the beginning 


_ of the calculations. Table 1 gives an abstract of the results discussed above. 


The aim of the present work is to obtain a molecular orbital (MO) description 
of the H,O* ion when unaffected by~charges of neighbouring atoms and ions. 


_ As is seen from the preceding, there are no definite values of angles and inter- 
atomic distances given as a starting point for such a calculation. The only ex- 


perimental values of these parameters have been obtained from H,0° ions in 


ionic crystals. These data must therefore deviate from those for a free ion and 
also must vary among themselves because they are taken from crystals of differ- 


ent complexes. In addition to this there are rather great uncertainties in the 


_ determination of the atomic positions of a given crystal, especially the positions 


of hydrogen atoms. For example, the standard deviation of determinations of 


. hydrogen positions with the X-ray method is about 0.07 A. However, from the 


results given in Table 1 it can be seen that the H-O-H angle preferably occurs 
in the interval 100°-120°. A tendency towards a more precise value is difficult 
to find. For a better understanding of the behaviour of the ion in its different 
associations we have chosen to study it assuming three different H-O-H angles: 
100°, 110°, and 120°. For a free H,0* ion it would, however, be justified to 


assume the equilibrium angle not smaller than in water, about 105°. 


Concerning the O-H distance, a study of Table 1 makes it seem convenient 
to choose a value in accordance with the statement by Richards and Smith [9]: 
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of the water molecule. Therefore» 
the O-H distance of water 0.9581 A 2 


~~. 


Il. Wave functions 


The system consists of one oxygen nucleus, O, three hydrogen nuclei Hy, Hy, 
and Hy, supposed to be arranged symmetrically, and ten electrons, whose con- 
figuration is to be determined (see Fig. 1). The symmetry group is C3, and for 
the flat molecule (H-O-H angle 120°) D3,. In the ground state the electron 
configuration is to be described by five doubly occupied molecular orbitals (MO’s) 
z ¢;. The total wave function is formed with the aid of a Slater determinant: 


a (1) J, (1) a (2) p; (2) a (3) 1 (3)... | 
B (1) fy (1) B (2), (2) B (3) 4, (3) «.. 


(1) 
a (1) $3 (1) a (2) ps (2) a (3) (3)... 


where 1, 2,3... denote the coordinates of the electrons. « and f# are functions 
representing the two possible spin orientations. Each MO is formed as a linear 
combination of atomic orbitals (LCAQ): 


pi = Ding Xx (2) 


where %, are atomic orbitals (AO’s) centered on the different atoms in the mole- 
cule. The atomic orbitals (see Fig. 1) used here are: 


Oxygen: 1s, 2s, 2px, 2py 2pz, abbreviated 0, s, x, Y, 2, . 
Hydrogen: 1s on each hydrogen atom denoted by hi, ho, hs. 


The 2s function has been made orthogonal to 1s by giving it a node according | 
to 2s=(1—S*) *(2s’—S-+1s), where 2s’ is nodeless. S is the overlap of 1s with 
2s’. Orbital exponents are given by Slater’s rules: oxygen 1s 7.7, 238’, 2pu, 2py, 


2pz 2.275, hydrogen 1s 1.0. The best coefficients a, can be determined by 
minimizing the energy 


—— 


E=(Y*HWdr (3) 


of the system described by ’. The hamiltonian H has the following form: 


10 2 Bi 
h 8 1 I iin. ah 
=— > |—A,+—+ =\+5 =, 
Py (* Age 2, Yan) 2 2 Yi (4) 
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rt problem is further treated according to the iteration method developed by 

oothaan [18]. To facilitate the solution of the secular equation which arises 
om the minimum condition it is better to express the MO’s as linear combinations 
yymmetry orbitals (SO’s) which are combinations of the AO’s earlier mentioned. 
$0O’s which form a basis for the irreducible representations of the symmetry 
p, can be obtained with the aid of Wigner’s projection operators [19]. The ; 
ansformation to the SO’s o; can be written: 


~ 


: 1 
1 05 == (hy — hs) 
= — (hy hy + hg) e_ V2 
V3 
@, ay; 
4 1 
03—8 G7 = —= (2h, — hy — he) 
i Cu kts e+ V6 
= 0,—%, 


_ where the symmetry types are written before each group. 
The MO’s are now obtained as linear combinations of the symmetry orbitals: 


e- 
di = (ny) = Seu o,- for 1=1...4 


ke 


i =(ne_) = Ci5 O05 + Cig 06 for i=5, 6 (6) 


di=(ne4)=Ciz07+Cig03 for 1=7, 8. 


* 


i = 
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The symmetrical ground state is represented by (ref. 18, p. 83): (1 a,)? (2a,)? 


(3a,)? (1e_)?(1e,)®. The orbitals belonging to e_ and e, symmetries are degenerate. 
Therefore the coefficients c,, and cs, will be equal to c,, and cz, respectively. 

In this way the secular determinant will reduce to one four-by-four and two 
two-by-two determinants. The latter two are equal because of the degeneracy 
mentioned above. 


II. Integrals 


As pointed out earlier, most of the integrals are obtained from the tables by 
Ellison and Shull [16]. The integrals. which they had to approximate in some 
way, they usually estimated by different methods; some integrals were also ex- 
panded by Ruedenberg’s method [20] to an almost exact value. A separate dis- 
cussion of molecular integrals has been given by Ellison [17]. However, according 
to McWeeny and Ohno [34], some of the integrals in Ellison and Shull’s paper 
are in error. With two exceptions the revised integrals deviate from the original 
ones only in the fourth decimal. The two integrals S(h’s) and (h’z’:ss) are 
revised from 0.4946 to 0.5086 and 0.2073 to 0.2193 respectively. We changed 
these two integrals in this way and repeated the calculations, after first having 
used those of Ellison and Shull. Coefficients and other quantities where somewhat 
changed thereby, but not sufficiently to affect the general conclusions drawn 
from them. The remaining integrals of three- and four-center type we have tried 
to compute with accuracy to the third decimal. This is described in more detail 
in the appendix. 


IV. Final coefficients and orbital energies 


The final coefficients making the iteration procedure self-consistent are listed 
in Table 2, together with one-electron energies in atomic units. 

As might be expected, the mixing between the oxygen 1s orbital o, and the 
others is small. However, it is not negligible. In their calculation on the water 
molecule Ellison and Shull [16] made, besides a complete calculation allowing 
such mixing, one other more approximate calculation where they did not allow 
mixing, letting the 1s electrons belong to the oxygen core. The results of these 
calculations differed by about 20% both in dissociation energy and dipole mo- 
ment. Discrepancies caused by neglect of 1s mixing in this case seem to be of 
about the same magnitude as the discrepancies in results obtained by conven- 
tional LCAO MO calculations, which are caused by the approximate nature of 
this method. 

The values of ¢, which give good approximations to the negatives of the 
corresponding ionization energies, show a definite difference between occupied and 
unoccupied orbitals in the entire investigated angle interval. The effect of con- 
figurational interaction would therefore be supposed small. 


V. Total energy and proton affinity of water 


_The total energy for the different angles appears in Table 3. The expression 
giving the energy as a function of the angle m (see Fig. 1) must be even in @ 
by symmetry arguments. The lowest polynomial fitting our computed values is 
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Table 2. Final coefficients c;, and orbitals energies ¢; for the three different values 
of the angle m (see Fig. 1). 


a) 100° 110° 120° 


SS Ce 


Cis — 0.0046 — 0.0045 — 0.0044 

Cra 1.0004 1.0004 1.0004 

4 C43 0.0178 0.0177 0.0176 
Cus 0.0027 0.0019 0 

Cox 0.0953 0.0848 0.0780 

Coa — 0.0233 — 0.0223 — 0.0217 

Bienupiea Cog 0.8867 0.9096 0.9286 
Cag 0.1947 0.14.00 0 
‘iad 0.2372 0.1694 0 
C59 — 0,0221 — 0.0159 0 
C33 — 0.4460 — 0.3169 0 
dee 0.8777 0.9397 1 

Css = Cr7 0.5397 0.5419 0.5495 

Csg =Crg 0.7269 0.7081 0.6880 

ad 1.0076 1.0458 1.0771 

Cas — 0.1013 — 0.1054 — 0.1087 

Gieotsio’ C4 — 0.9236 — 0.9739 — 1.0209 
Cas — 0.5264 — 0.3748 0 

Ces = Cer 1.3683 1.3341 1.3091 

Ces = Cas — 0.8826 — 0.9140 — 0.9448 

& — 21.0608 — 21.0468 — 21.0410 

Ovugiea bs — 1.8440 = 1.8171 — 1.7944 

&5 — 0.9718 — 0.9516 — 0.9392 

&,= £5 — 1.0986 — 1.1156 — 1.1356 

Vaeccapica & 0.0090 0.0651 0.1108 

£5 =&s 0.1025 0.0947 0.0675 


Table 3. Total molecular energy for H,0* and proton affinity of water. 


@ 100° Lio? 120° Exp. 
Total energy in a.u. — 76.0836 — 76.1407 — 76.1820 — 76.725 
Proton affinity of 
water in kcal/mole 202 238 264 169 [27] 


then of fourth degree. This function has one extremum in the interval: a mini- 
mum at m»=0, which means w=120°. The experimental data mentioned in the 
introduction indicate the possibility of obtaining a minimum point just below 
120° and a flat maximum at 120° (see discussion, part I). 

The value of the proton affinity, which is the heat evolved in the reaction 
H,0 + H*—H,0%, is given as the difference between the total energy of the water 
molecule and the H,0* ion. For a theoretical calculation the most consistent way 
to choose the water energy is in this case to take the value calculated by the 
same method and with the same basic orbitals as for H,0". The energy at 
the experimental equilibrium angle is according to McWeeny and Ohno [34] 
— 2061.0 e.v.= — 75.761 a.u. The difference between this and the total energies 
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of H,O* are 202, 238, and 264 kcal/mole, for the angles 100°, 110°, and 120° 
respectively (Table 3). 

The experimental value appears to lie between 150 and 200 kcal/mole. Sherman 
[24] calculated it with the thermochemical egpaee 


Pu,0 = Pyu, + Qn, — Qu,0 + Qu,oc10, — Qnx,c10, 


and obtained the result 182 kcal/mole. The value of Sherman corrected to 188 
kcal/mole [25] has been the value generally accepted. Later on Vetchinkin, 
Pshenichnov, and Sokolov [26] pointed out that if a correction is made account- 
ing for the energies of the hydrogen bonds related to the hydrogen atoms in 
H,O* instead of making rougher models of the ion in the calculations of some ~ 
lattice energies, the energy will be lowered to about 170 kcal/mole. This is in very — 
good agreement with the value obtained by Tal’rose and Frankevich [27]—169 + 2 
kcal/mole—with the mass spectroscopic method (see also ref. 28). 

The affinity obtained from our theoretical investigation may be said to be 
satisfactory, especially since the total affinity is only 0.55% of each of the two 
quantities whose difference it is. 


VI. Equivalent orbitals 


The MO’s in (6) have been transformed to equivalent orbitals (EO’s) [29]. In 
this case the EO’s will consist of three equivalent bonding orbitals %»1, X52, and 
%v3, one lone-pair electron orbital %, and one oxygen inner shell orbital. The 
last is supposed to constitute the molecular orbital ¢,. The four remaining MO’s 
do, b3, $5, and ¢, then undergo a unitary transformation to the four equivalent 
orbitals %o1 ... Y»3 and X,. This transformation is not uniquely determined by the 
requirement that the %’s are orthogonal and equivalent. We can put one more 
condition upon them. The equivalent orbitals for two different forms of this con- 
dition have been evaluated: 


(i) The lone-pair electron orbitals X, shall contain only oxygen atomic orbitals. 
(ii) The oxygen hybrid part of the bonding orbitals shall be directed along the 
O-H. direction. 


As can be seen from Table 4, the two conditions do not differ very much 
from one another. In case (i) the bonding hybrid is put further away from the 
lone-pair hybrid than in case (ii). This inward bending of the bonds seems natu- 
ral in view of the repulsion from the lone pair. The contribution from the hy- 
drogen orbitals not belonging to the bond, b;2, is weakly outward-bending, but 
cannot affect the charge distribution much. 

In case (ii) we have forced the O-H direction upon the oxygen part of the bond- 
ing orbitals. The coefficient 6, which in case (i), @=100° and 110°, was small 
and of little importance for the charge density, has in case (ii) values indicating 
strong sp hybridization. Even the lone-pair orbital is in this case strongly hy- 
bridized. If we want bonding and lone-pair orbitals to have the property of 
pointing as much as possible away from each other, then picture (ii) is the most 
satisfying. This property is even more pronounced because the negative contri- 
bution J, causes the lone pairs to be still more removed from the O-H region. 
Such mixing in lone-pair orbitals can be of importance in keeping them away 
from regions “reserved” for bonding electrons. 
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Table 4. Equivalent orbitals of H,O* for different values of the H-O-H angle. 


(i) no hydrogen contribution to x 1 (ii) O part of x,, along O-H, line 
@=100° @=110° @=120° @=100° @=110° @=120° 
re i ne ee EY 


%p1 and x, 


0.5935 0.5781 0.5617 b, 0.5935 0.5781 0.5617 
—0.5121 —0.5214 0 b, — 0.3130 0.1984 0 

— 0.0479 0.0715 0.5361 db, 0.3585 0.4720 0.5361 

—0.0169  —0,0140  -0.0125 b, —0.0176 -0.0146  -0.0125 

0.4450 0.4244 0.3923 Dat 0.4216 0.4013 0.3923 

Sued? PF O7175 "0.86872" | °b)., be) 0.1181 “= 0.1406" 120.1572 
0.1467 0.2954 1 ke 0.7172 0.8858 1 
0.9891 0.9553 0 l, 0.7743 0.5095 0 
—0.0134  — 0.0128 0 ie 0.0104 0.0106 0 
0 0 0 bs —0.1016  —0.0846 0 


The coefficients b,, b,, ete. are the coefficients of the AO’s x, z, etc. in the equi- 
valent bonding orbital %,1. The two other bonding orbitals are obtained through 
rotation of %); +120°. The lone- -pair orbital coefficient 1, is the coefficient of each 
of the AO’s h,, h,, and h, in EO %;. The notations (i) and (ii) refer to the two 
cases discussed in the text 


VII. Population analysis 


Population analysis has been performed according to the scheme of Mulliken 
[30]. The results can be seen in Table 5. The notations in our special case may, 
however, be given some extra explanation: 


Partial gross population in MO ¢; and AO h,: N (i; h,). ‘ 
Partial gross population in MO ¢; and all hydrogen atoms: N (i; 3H) =e 
N (1; h,). 


N(t; 3H 
Total gross population on one hydrogen atom: N aay) 


, where the 


summation is to be performed over all occupied MO’s. 
Net charge of one H atom: @ (H)=1—WN(H). . 
Partial overlap populations between SO o, and all hydrogen AO’s in MO 4;: 


- 3H 4.41 Ci Sip for r= ey a 

be n= ; 
tes A Oop Cp Ong = 4 Cy7 Or, Sra LOK = 5; 7. 
Subtotal overlap population between SO o, and all hydrogen AO’s: 


n(3H, r)=>, n(t, 3H, 1). 


Total overlap population between oxygen and all hydrogen atoms in MO 4: 


n(i, 3H, O) => n(i, 3H, r) 


es 


: awe 
where the summation runs over the oxygen SO’s. 
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Table 5. Population analysis in the MO’s and SO’s of (6). Notations: see above 
and ref. 30. 


@ 100° HIQ® 120° 


ee ee ae ing ee ee ee oe eee ee 
N (1; 3H) — 0.0010 — 0.0010 — 0.0010 
N (2; 3H) 0.1915 0.1648 0.1472 
N (3; 3H) 0.1311 0.0653 0 
N (5; 3H)=N (7; 3H) 0.6477 0.6879 0.7313 
N (a,; H) 0.1072 0.0764 0.0487 
N (e; H) 0.4318 0.4586 0.4875 
N(H) 0.5390 0.5350 0.5363 
N (1; 2) 2.0005 2.0005 2.0005 
N (2; 2) 0.0006 0.0006 0.0006 
N (3; 2) — 0.0001 — 0.0001 0 
N (2) 2.0010 2.0011 2.0011 
N (1; 3) 0.0005 0.0005 0.0005 
N (2; 3) 1.7216 1.7908 1.8523 
N (3; 3) 0.2114 0.1063 0 
N (3) 1.9335 1.8975 1.8527 
N (1; 4) 0.0000 0.0000 0 
N (2; 4) 0.0863 0.0439 0 
N (3; 4) 1.6576 1.8285 2 
N (4) 1.7439 1.8723 2 
N (5; 6) =N (7; 8) 1.3523 1.3121 1.2687 
n(1, 3H, 2) — 0.0019 — 0.0019 — 0.0018 
n (2, 3H, 2) — 0.0009 — 0.0008 — 0.0007 
n(3, 3H, 2) — 0.0022 — 0.0011 0 
n (3H, 2) — 0.0051 — 0.0038 — 0.0026 
n(1, 3H, 3) — 0.0003 — 0.0003 — 0.0003 
n(2, 3H, 3) 0.2979 0.2718 0.2551 
n(3, 3H, 3) — 0.3727 — 0.1892 0 
n (3H, 3) — 0.0751 0.0823 0.2548 
n (1, 3H, 4) 0.0000 0.0000 0 
n(2, 3H, 4) 0.0209 0.0093 0 
n(3, 3H, 4) 0.2340 0.1245 0 
n (3H, 4) 0.2549 0.1338 0 
n(1, 3H, O) — 0.0022 — 0.0022 — 0.0021 
n (2, 3H, O) 0.3179 0.2803 0.2544 
n(3, 3H, O) — 0.1409 — 0.0658 0 
n (3H, O) 1.3575 1.4495 1.5408 
n(a, 3H, O) 0:1747 0.2123 0.2522 
n(5, 3H, 6)=n(7, 3H, 8) 0.5914 0.6186 0.6443 
Q (H) 0.4610 0.4650 0.4637 


Total overlap population between oxygen and all hydrogen atoms: 


n(3H, 0)=> n(i, 3H, 0). 


As might be expected the net charge Q(H) of the hydrogen atoms is strongly 
positive, about +0.46 for all the angles considered. The molecular orbitals of ay 
symmetry, ¢,, ¢,, and ¢, have small contributions from the hydrogens. The total 

3 - sie 


hydrogen population from them, N (a,,3H)=> N (i, 3H), varies between 0.15 and 
¢ : = 
about 0.32, while that from the e orbitals ¢, and ¢; lies between 1.29 and 1.46. 
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The strong oxygen-hydrogen mixing in the e orbitals also make them the strongest 


bonding orbitals. The total overlap population between oxygen and all hydrogen 
atoms in each of the MO’s ¢; and ¢,, n(5, 3H, 6), is twice the corresponding 


population, n(2,3H, 0), in the bonding a, orbital. The orbital $3 is antibonding, 
but the overlap population n(3,3H, O) decreases to zero by symmetry as the 


_ H-O-H angle increases to 120°. The relatively weak bonding of the orbital gd, may 


- 


be explained by the fact that the molecule for all the treated angles is rather 
flat, which diminishes the possibility of the only oxygen 2p function of a, sym- 
metry, o4, joining the 2s function to form a strong 2s—-2p hybrid that could make 


_ & great contribution to the bond. Therefore, the bonding is predominantly of the 
_ weaker 2s—ls type. 


VIII. Electric moment 


The center of gravity of the charges is in both the non-flat cases situated some- 
what further from the oxygen atom than the center of mass (see Table 6.) This 


_ electric moment influences the rotation and vibration spectra of the molecule ion. 
_ Any experimental value of the moment has not been found (see discussion). 


Table 6. The electric moment of the H,O* ion as a function of the H-O-H angle. 
In Debye units. 


@ | 100° KOS 120° 


Electric moment | 2.96 Paella) 0 


IX. Discussion 


If an experimental value of the electric moment of the H,O~ ion could be ob- 
tained, it is probable that, in connection with our theoretical study of its depen- 
dence on the H-O-H angle, it would give a rather satisfying answer to the 
question of the real angle. Molecular orbital descriptions usually give passable 
values of the dipole moments. In addition, the moment in our case varies very 
rapidly with the angle in the interval 110°-120°, which would make it easier to 
fix the real value of the angle. In their study of the H,0* ion in crystals of 
hydrogen halides Ferriso and Hornig [7] found a possible indication of the elec- 
tric moment. They could not find the band corresponding to the torsional vi- 
brations of the H,O* ions (librations) perpendicular to the symmetry axis. They 
concluded: ‘‘If the bands at 770 cm‘ represent this frequency, the dipole moment 
of the OH} ion must be very small, at least ten times smaller than in ammonia. 
If this is true, the moment would be smaller than 0.15 D, which, according to 
Table 6, would indicate an H-O-H angle only negligibly different from 120. 

This equilibrium angle is also the result of our study of the angular depend- . 
ence of the total energy. Since, however, the variable part of the total energy 
is very small compared with the total energy itself, we cannot assign this result 
any great significance, especially since we have let the O-H distance remain fixed 
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under the angle variation. In ammonia, NH,, with which H,O* is isoelectronic, 
the equilibrium H-N-H angle is 107°. In H,O*, where the O-H distance is about 
the same as the N-H distance in NH,, one can expect the angle to be definitely 
greater because of the excess charge of the ion, which appears mainly on the 
hydrogen atoms, thus causing a considerable extra repulsion between them. Con- 
cerning the H,O* ion in the experimentally studied crystals, it is probable that 
the rather different angle values which have been proposed correspond to dif- 
ferent values even in reality. The large positive charges of the hydrogen atoms 
may be a very effective cause for strong bending forces from the frequently 
ionized neighbouring atoms in the ionic crystals studied. In heavy ice, with a 
much smaller positive charge on the deutorium atoms, the D-O—D angle is within — 
a half degree of the tetrahedral angle and thus seems to be entirely subordinate 
to the surrounding crystal structure, which has tetrahedral angles [31, 32]. It is 
possible that in a future paper we shall be able to obtain more detailed results 
of the interaction of the hydronium ion with its lattice neighbours. We are es- 
pecially interested in the properties of the hydrogen bonds where hydronium ions 
take part. Besides the appearance in ionic crystals like those mentioned earlier, 
the ion is thought to appear momentarily in hydrogen-bonded systems such as 
ice, where they are intermediate states in the proton transport process, which is 
thought to consist partly of proton jumps from one position in a hydrogen bond 
to the other [2, 3, 33] 
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APPENDIX 


Three-center integrals 


The coulomb interaction one-electron integral (H,:h,h,) and the two-electron 
integrals (hyhy:hgh3) and (hyh,:h,h3) have been computed in the following way. 


(H, : Ay hs) 


This integral has been obtained with four-decimal accuracy by an expansion 
given by Coulson [21]. It has then been used as the basis for the estimation of 
other integrals. We have also used other methods in order to compare them and to 
test their reliabilities in giving values for other integrals; theMulliken approximation 
[22], where the charge distribution y,y, is replaced by 2 Spy. (YrY1t YoY), in 
this case will clearly give a value of the integral which is too low, because it con- 
centrates charge around the two nuclei H, and H, instead of having constant 
density on the line H,-H;, which is the case here. Along this line the density 
has its maximum and will decrease rather rapidly beyond it. A 1s, 1s, distribution 
where a and 0b are similiar nuclei therefore may be fairly well represented by a 
homogeneously charged rod ab with total charge Ss, 1.,. This approximation gives 
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in our case a somewhat low value (Table 7). The integral has also been ap- 
‘proximated with a function of the angle 6= \ H,H,H, in such a way that Sr 
comes in as a linear factor: me 


- 


om 


¥ be dni 1p) =a 8 (p) +0. (7) 


For @=0 the integral reduces to a two-center one and for p= it can be cal- 
‘culated exactly according to Barker and Eyring [35]. Then a and 5b are fixed. 


Table 7. The integral (H,:h,hs) in different approximations. w= 110°. 


Method (H, : Aghs) 
Exact 0.1299 
Mulliken approximation 0.1185 
I=aS+b 0.1245 
Rod approximation 0.1254 


(hy hy: hy hs) 


The Mulliken approximation has been assumed to give the same percentage 
error in this integral as in (H,:h,hz,). 


(hy hy: hy hg) 


The Mulliken approximation of one of the distributions gives two- and three-cen- 

- ter integrals of types earlier treated. The same approximation of both the distribu- 

tions gives only two-center integrals. The results of the two methods for w= 110° 

are 0.0516 and 0.0503 respectively. Here it is more difficult to see if the ap- 

proximate values, which agree fairly well with each other, lie below or above the 
correct value. The value 0.0516 is probably the most reliable of them. 


Four-center integrals 
(h, 0: hghs) 

The Mulliken approximation of h,0, which places equal charge on H, and O, can- 
not be suitable for this very unsymmetrical distribution. The refinement of it, the 
Léwdin approximation [23], consisting in putting py, py .=1,yY1Y1+42Y2Y2, where 
4, and A, are determined to give the correct total charge and charge moment, 
is more suitable here. The result is in good agreement with that obtained by the 
semi point-charge approximation (SPCA) where the h,o distribution is approxi- 
mated by a point charge of magnitude S;,,. at the O nucleus (see Table 8). 


Table 8. Approximate values of the integral (h,o:h,h3). o=110°. 


Method (hyo 2 Ayhs) 
Mulliken approximation 0.0111 
Léwdin approximation 0.0127 
SPCA 0.0129 


[i se 


159 


R. GRAHN, A theoretical study of the H,O* ion 


(hy 8: hy hs) . 

The Mulliken approximation is supposed to give values of this and the integrals 
earlier studied (h,0:hgh3) and (hyh,:hgh3), which deviate from the exact ones by 
a factor which depends linearly upon the distance from the centers of gravity of 
the two charge distributions in question. This corrected Mulliken approximation 
gives the result 0.0841 for w= 110°. The Léwdin approximation of h,s may be — 
supposed to give a value which is too high because in it charge is moved from 
the region between the O and H, nuclei to the O nucleus, which effect pre- 
dominates over the contrary effect caused by the corresponding movement of © 
charge to the H, nucleus (see Table 9). 


. 


Table 9. Approximate values of the integral (hys:h,h3). o= 110°. 


Method (his : hghg) 
Mulliken approximation of hah, 0.0737 
Mulliken approximation of h,s 0.0809 
Corrected Mulliken approximation 0.0841 
Loéwdin approximation 0.0884 


(hy 2: hy hs) 


The approximation function 
I(p)=6 sin pS(¢), (8) 


where @ is defined by Fig. 1 and S(p)=Sh),,,, has the correct value at p=0 and 
correct derivative (=0) at g=2/2. We can choose b=IJ(z/2) to give J the cor- 
rect value at y=2/2. Its derivative at g=0 is not zero, which would be the 
case with a function of the type earlier used in (7). 

For p=110° the value obtained is —0.0300, in very good agreement with the 
value — 0.0302 obtained by the Mulliken approximation and corrected with the same 
factor as for (h,s:hgh3) (see Table 10). 


Table 10. Approximate values of the integral (h,z:h,h3). o =110°. 


Method | (A122 Aghs) 
Mulliken approximation | — 0.0265 
Corrected Mulliken approximation — 0.0302 
b sin PSh,n, — 0.0300 


(hy @: hy hs) 


The Mulliken approximation with the same correcting factor as for (h,s:h,h,) 
and (h,z:h ahs) has been used. aig 

Values of the above-mentioned integrals used in our calculations for the dif- 
ferent @ values are collected in Table 11. 
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Table 11. Values of some three- and four-center integrals. (w defined in Fig. 1.) 


a 
} 


- 


Integral \ 100° 110° 120° 
(Hy :hghs) | 0.1535 0.1299 0.1121 
(hyh, : Rg hs) 0.1472 0.1257 0.1092 
(hy hy : hy hg) 0.0663 0.0516 0.0411 
(hyo : hg hg) 0.0141 0.0129 0.0119 
(ys : hg hg) 0.0948 0.0841 0.0752 
(hyz:Rghs) | —0.0480 —0.0300 0 
(hya : hgh) 0.0254 0.0258 0.0264 
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